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Abstract. A populations of stellar mass black hole binaries may exist in 
globular clusters. The dynamics of globular cluster evolution imply that there 
may be at most one black hole binary is a globular cluster. The population of 
binaries are expected to have orbital periods greater than a few hours and to 
have a thermal distribution of eccentricities. In the LISA band, the gravitational 
wave signal from these binaries will consist of several of the higher harmonics of 
the orbital frequency. A Monte Carlo simulation of the galactic globular cluster 
system indicates that LISA will detect binaries in 10% of the clusters with an 
angular resolution sufficient to identify the host cluster of the binary. 
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1. Introduction 

Globular clusters may be the first isolated stellar systems to have their population of 
relativistic binaries studied with gravitational radiation. The galactic globular cluster 
system consists of about 150 clusters scattered throughout the halo of the galaxy 
At an age of nearly 13 Gyr, they consist of old population II stars and stellar remnants. 
The dynamical evolution of a globular cluster will tend to concentrate the more massive 
stars and binaries into the core of the cluster. This will enhance the likelihood of 
encounters between binaries and single stars, which tends to increase the number of 
relativistic binaries in globular clusters. These binaries will be potential sources for 
LISA §. 

Although there is no direct observational evidence for black hole binaries in 
globular clusters, there are two theoretical scenarios which are expected to produce 
black holes in globular clusters. Miller and Hamilton describe a mechanism whereby 
intermediate mass black holes will be produced in globular clusters with high central 
densities. These black holes will have masses ~ 1000 Mq and may reside in the cores 
of up to 20 globular clusters. As compact objects in the core of these globular clusters 
pass close to the intermediate mass black hole, they may emit gravitational radiation 
within the frequency band of LISA. This radiation will be in the form of bursts. There 
is currently evidence of an intermediate mass black hole in M15 Q. 

The scenario of Portegies Zwart and McMillan allows for the production of 
numerous stellar mass black hole binaries in the cores of globular clusters. These are 
then ejected from the cluster through interactions with other black holes and binaries 
in the cluster. The last black hole binaries in the cluster will have no other black holes 
to eject it. Thus, the last dynamical encounter can have three possible outcomes: the 
binary and the encountering black hole can all be ejected, the binary can be ejected 
and the encountering black hole can be retained, the binary can be retained and 
the encountering black hole can be ejected. Although these three outcomes are not 
equally likely, the probability of retention of a black hole binary can be estimated to 
be between 30% and 50%. Consequently, between 30% and 50% of globular clusters 
can be expected to have one remaining stellar mass black hole binary. This highly 
eccentric binaries of this population are viable sources for LISA. 

2. Stellar Mass Black Holes 

Stellar mass black holes with M ~ 10 Mq are formed in globular clusters within the 
first 10* yr of its history. These objects will segregate to the core of the cluster within 
the first few relaxation times. Barring any massive seed black holes, these objects 
will be the most massive objects in the core. Consequently, they will effectively form 
their own cluster and begin to interact solely with themselves in the presence of the 
background potential of the cluster. Encounters with the other stars in the cluster will 
provide little more than slight perturbations to the subsequent dynamical evolution of 
the black holes. This evolution will lead to the ejection of nearly all of the black holes 
in the cluster. The ejected black holes were studied as possible sources for LIGO by 
Portegies Zwart and McMillan Q. However, the dynamical evolution may also leave 
a few longer period black hole binaries in the cores of globular clusters. These are the 
potential sources for LISA. 

The process by which dynamical evolution leads to the ejection of black holes 
from globular clusters also leads to the hardening or shrinking of the binary systems. 
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The outcome of a gravitational encounter between a binary and a field star in a 
cluster depends upon the binding energy of the binary and the masses of all three 
bodies involved in the encounter. If the binding energy of the binary is greater than 
the kinetic energy of the field star, then the common result is a transfer of binding 
energy to the kinetic energy of the three stars involved in the encounter. The least 
massive of the three stars is ejected from the encounter and the remaining two form 
a tighter binary. The binding energy of the system generally increases by about 20% 
in an encounter and roughly 1/3 of the liberated gravitational energy goes into recoil 
energy of the binary. The process results in the gradual ejection of both binary and 
single black holes from the globular cluster as they gain sufficient velocity to escape 
the potential of the cluster. 

This process is expected to take roughly 2 Gyrs. At the end of this time, the 
most likely result is three black holes remaining in the form of a binary and a single 
black hole. The last binary/single interaction can result in three possible outcomes: 
all three black holes can receive sufficient recoil to escape from the cluster, a binary 
can remain and eject a single black hole, or the single black hole can remain and eject 
the binary. The last outcome is highly unlikely as it would require the single black 
hole to be more massive than the binary. The more likely outcome would then be 
the exchange of the massive black hole with the least massive of the original binary 
components, resulting in the ejection of the least massive black hole and retention 
of the newly created massive binary. Thus, we can expect that roughly half of the 
globular clusters to contain a retained black hole binary. 

The binaries which remain in the globular clusters will have a thermal distribution 
of eccentricities so that P{e) = 2e. We can determine the minimum orbital period 
expected for a retained binary by looking at the maximum orbital period of an ejected 
binary. Following Portegies Zwart and McMillan, a binary will be ejected if 

El. > ^iGWo^kT (1) 

[m) 

where Wq = {m]^\^)\(f>()\/kT is the dimensionless central potential of the cluster and 
(3/2)fcT is the average kinetic energy of a cluster star. The average cluster star mass 
is (m) and mbh is the mass of each black hole in the binary. If the cluster is in virial 
equilibrium, then kT = GM^ /GNrvh- where M is the total mass of the cluster and 
Tvir is the virial radius. Consequently, the binding energy of a binary with semi-major 
axis a can be related to the descriptive parameters of the cluster by: 

E^^SNl^Y'^kT. (2) 
Combining Equations Q and || with Kepler's third law gives 

for the minimum orbital period of any remaining black hole binary in a globular 
cluster, where we have assumed that N{m) — M. 

We can estimate the range of periods for the remaining black hole binaries by 
inserting standard values for cluster parameters into Equation ^. Using estimates from 
Takahashi and Portegies Zwart Q of the parameters of zero-age globular clusters, we 
use logrvir = ± 0.3pc, logM = 6.0 ± 0.5Mq, and Wq = 7.5 ± 2.5 to determine this 
range. We find that the Porb ~ 900 — 120, 000s with the average value of Porb — 4, 200s 
for black hole masses of mbh = lOM©. Although this range is for the minimum orbital 
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period of the remaining black hole binaries, it is not unreasonable to assume that the 
actual orbital periods will be close to the minimum value. 

3. Signal at LISA 

In a coordinate system centered on an eccentric binary with orbital frequency /, the 
gravitational wave at a distant detector located at angular position i) and tp is given 
by the two polarizations |7| : 

/ix [^^^y^ (27m/t) cos2(^ - h''^]!y cos (27m/t) sin 2^ (4) 

= ^ E + cos' ^) sin (2^n/<) sin 2^ 

^ n— 1 

+ h''^}y COS (27m/t) cos 2iy9 
- (1 - cos^ cos (27rn/i) , (5) 

where the metric components h''^ and /i^'^j, are: 

h<-^J - /lo?! (1 - e^) [ J„_2 (r^e) - 2 J„ (ne) + J„+2 (ne)] (6) 



r 2 

ft-i'L^j, = 2honl^Jn-2ine) - 2eJ„_i(ne) + -J„{ne) 

+ 2eJ„+i(ne) - J„+2("e)| (7) 

= -4/ioJ„(ne). (8) 
The common amplitude factor ho is: 

ho = ^(2^/)2/3a1^/3 (9) 

with the "chirp mass" A^^/^ = MiA^iMi + Mj)"^/^. 

In the frequency domain, such a signal will appear in LISA as a series of 
monochromatic signals, each of which is modulated by the motion of the detector 
about the Sun. We can describe the signal received by LISA in terms of several 
defined unit vectors. Let the orientation of the arms of the constellation of the LISA 
spacecraft be given by the three vectors, £i, £2, is- The direction from the sun to 
the source is n. The orientation of the binary orbit defined by L (which points along 
the angular momentum vector), a (which points along the semi-major axis toward 
periastron), and b = L x d. The gravitational wave propagates along —n, and we 
define the polarization axes p = (^n x Lj j\ri x X| and q = —nxp. The angles used in 

Equations ^ and ^are then defined by cos d — —fi L and tan ip = (^n ■ bj / (n ■ d). The 
sensitivity of one interferometer of LISA to both polarizations can then be described 
by F+ a,nd F"", with: 

F+ = {paPb - qaqbKnel ~ e^24) (10) 

F'^ = {paqb-qaPbWA-r2i'2)- (11) 
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The signal received at LISA is then: 
hit) 



OO 

E 

n=l 



^/Al + Bl cos [27rn/i + 0„p(t) + 0„d(O] 



(12) 
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(13) 

(14) 
(15) 
(16) 

The angular position of the source in heliocentric coordinates is thetas and 0s- The 
center-of-mass trajectory of LISA is 0(i) = 2TTt/T with T = 1 yr and R = 1 AU. A 
second signal can be synthesized from a linear combination of the responses from all 
three arms of LISA. This signal is equivalent to the signal from a detector rotated by 
45° in the plane of the initial detector 

We have performed a study of the gravitational wave signal from a Monte Carlo 
simulation of the population of black hole binaries in globular clusters. We have chosen 
a pessimistic retention probability of 1/3, so that there are 43 black hole binaries in 
the galactic globular cluster system. The black hole masses were chosen randomly 
from the range 6 — 18Mq. The eccentricities were chosen from a thermal distribution 
(P(e) = 2e, and the orbital period was selected from a flat distribution in the range 
2, 000 — 100, 000s. The angular resolution of the gravitational wave signals from these 
black holes was determined by looking at the first 10 harmonics above 1 mHz. These 
black hole binaries were assigned random globular cluster hosts, and the signal-to- 
noise was then calculated using the expected instrument noise and transfer function 
for LISA calculated by Larson and Hiscock ||]. The results show that fourteen of 
these had signal-to-noise greater than 10. Figure |^ shows the eccentricities and orbital 
periods for all 43 black hole binaries in the simulation. 

The Fisher matrix formalism to find the expected variances for the 
parameterization of the signal from eccentric binaries allows us to find the angular 
resolution for the binaries in the simulation. Using only the first five harmonics above 
1 mHz, we find that the fourteen binaries with signal-to- noise above 10 also have 
sufficient angular resolution to identify them with their host clusters. The location 
and angular resolution for these 14 binaries are shown with the galactic globular cluster 
system in Figure 0. 



4. Conclusions 



The galactic globular cluster system may contain numerous black hole binaries. The 
evolution of stars in the early life of a globular cluster nearly guarantee that there 
will be several stellar mass black holes in each globular cluster. Through dynamical 
interactions and mass segregation, these black holes are expected to form a separate 
population within the cores of globular clusters. Further dynamical interactions 
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Figure 1. Signal-to- noise values for all black hole binaries in the simulation. The 
data points are plotted for orbital period (in seconds) and eccentricity. 




Figure 2. Galactic globular cluster system in heliocentric coordinates. Black hole 
binaries are shown as gray squares, while the globular clusters are black triangles. 
The angular resolution for black hole binaries with error boxes small enough to 
allow for identification of the host cluster is shown. 
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between the black holes will generate a population of black hole binaries which will 
subsequently be ejected from the globular cluster through recoil. However, it is 
expected that the final three black holes in the cluster will either eject each other, 
or result in a remaining binary system. During the interaction, the probability of 
retention can be estimated to be somewhere between 1/3 and 1/2. These retained 
binaries will have orbital frequencies below the gravitational wave frequency band of 
LISA. However, they are also expected to have large eccentricities. The luminosity 
in gravitational radiation from these binaries will be spread out among many of the 
higher harmonics, and consequently into LISA frequency band. We have performed a 
Monte Carlo simulation showing that ~ 10% of the globular clusters may contain a 
black hole binary with sufficient signal strength in the higher harmonics to allow for 
the angular position of the binary to be resolved to within the host cluster. Because 
the period distribution and likelihood of ejection for the black hole binaries are related 
to the structure and past history of the cluster, the identification and study of these 
black hole binaries may provide information about the past structure and history of 
their host clusters. 
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